From Milliseconds to Minutes: Melittin Self-Assembly from Concerted Non-Equilibrium Pressure-Jump and Equilibrium Relaxation Nuclear Magnetic Resonance.

Non-equilibrium kinetics techniques like pressure-jump nuclear magnetic resonance (NMR) are powerful in tracking changes in oligomeric populations and are not limited by relaxation rates for the time scales of exchange that can be probed. However, these techniques are less sensitive to minor, transient populations than are Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments. We integrated non-equilibrium pressure-jump and equilibrium CPMG relaxation dispersion data to fully map the kinetic landscape of melittin tetramerization. While monomeric peptides weakly form dimers (Kd,D/M ≈ 26 mM) whose population never exceeds 1.6% at 288 K, dimers associate tightly to form stable tetrameric species (Kd,T/D ≈ 740 nM). Exchange between the monomer and dimer, along with exchange between the dimer and tetramer, occurs on the millisecond time scale. The NMR approach developed herein can be readily applied to studying the folding and misfolding of a wide range of oligomeric assemblies.

Experimental NOE, Chemical Shift, and Proline Isomerization Data Provide Detailed Insights into Amelotin Oligomerization.

Amelotin is an intrinsically disordered protein (IDP) rich in Pro residues and is involved in hydroxyapatite mineralization. It rapidly oligomerizes under physiological conditions of pH and pressure but reverts to its monomeric IDP state at elevated pressure. We identified a 105-residue segment of the protein that becomes ordered upon oligomerization, and we used pressure-jump NMR spectroscopy to measure long-range NOE contacts that exist exclusively in the oligomeric NMR-invisible state. The kinetics of oligomerization and dissociation were probed at the residue-specific level, revealing that the oligomerization process is initiated in the C-terminal half of the segment. Using pressure-jump NMR, the degree of order in the oligomer at the sites of Pro residues was probed by monitoring changes in cis/trans equilibria relative to the IDP state after long-term equilibration under oligomerizing conditions. Whereas most Pro residues revert to trans in the oligomeric state, Pro-49 favors a cis configuration and three Pro residues retain an unchanged cis fraction, pointing to their local lack of order in the oligomeric state. NOE contacts and secondary 13C chemical shifts in the oligomeric state indicate the presence of an 11-residue α-helix, preceded by a small intramolecular antiparallel ß-sheet, with slower formation of long-range intermolecular interactions to N-terminal residues. Although none of the models generated by AlphaFold2 for the amelotin monomer was consistent with experimental data, subunits of a hexamer generated by AlphaFold-Multimer satisfied intramolecular NOE and chemical shift data and may provide a starting point for developing atomic models for the oligomeric state.

Synergism between x-ray crystallography and NMR residual dipolar couplings in characterizing protein dynamics.

The important role of structural dynamics in protein function is widely recognized. Thermal or B-factors and their anisotropy, seen in x-ray analysis of protein structures, report on the presence of atomic coordinate heterogeneity that can be attributed to motion. However, their quantitative evaluation in terms of protein dynamics by x-ray ensemble refinement remains challenging. NMR spectroscopy provides quantitative information on the amplitudes and time scales of motional processes. Unfortunately, with a few exceptions, the NMR data do not provide direct insights into the atomic details of dynamic trajectories. Residual dipolar couplings, measured by solution NMR, are very precise parameters reporting on the time-averaged bond-vector orientations and may offer the opportunity to derive correctly weighted dynamic ensembles of structures for cases where multiple high-resolution x-ray structures are available. Applications to the SARS-CoV-2 main protease, Mpro, and ubiquitin highlight this complementarity of NMR and crystallography for quantitative assessment of internal motions.

Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings.

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, Mpro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic Rfree, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55-2.19 Å) Mpro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member "super ensemble" averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB Mpro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.

SARS-CoV-2 omicron variants harbor spike protein mutations responsible for their attenuated fusogenic phenotype.

Since the emergence of the Omicron variants at the end of 2021, they quickly became the dominant variants globally. The Omicron variants may be more easily transmitted compared to the earlier Wuhan and the other variants. In this study, we aimed to elucidate mechanisms of the altered infectivity associated with the Omicron variants. We systemically evaluated mutations located in the S2 sequence of spike and identified mutations that are responsible for altered viral fusion. We demonstrated that mutations near the S1/S2 cleavage site decrease S1/S2 cleavage, resulting in reduced fusogenicity. Mutations in the HR1 and other S2 sequences also affect cell-cell fusion. Based on nuclear magnetic resonance (NMR) studies and in silico modeling, these mutations affect fusogenicity possibly at multiple steps of the viral fusion. Our findings reveal that the Omicron variants have accumulated mutations that contribute to reduced syncytial formation and hence an attenuated pathogenicity.

Recombinant Expression and Chemical Amidation of Isotopically Labeled Native Melittin.

Post-translational modifications are ubiquitous in the eukaryotic proteome. However, these modifications are rarely incorporated in NMR studies of eukaryotic proteins, which are typically produced through recombinant expression in E. coli. Melittin is the primary peptide in honey bee venom. Its native C-terminal amide significantly affects its equilibrium structure and dynamics in solution and is thus a prerequisite for studying its native structure and function. Here, we present a method for producing triply isotopically labeled (2H, 13C, and 15N) native melittin through recombinant expression followed by chemical amidation. We then show that structural models produced with AlphaFold-Multimer are in even better agreement with experimental residual dipolar couplings than the 2.0 Å resolution X-ray crystal structure for residues G3-K23.

NMR Observation of Sulfhydryl Signals in SARS-CoV-2 Main Protease Aids Structural Studies.

The 68-kDa homodimeric 3C-like protease of SARS-CoV-2, Mpro (3CLpro /Nsp5), is a key antiviral drug target. NMR spectroscopy of this large system proved challenging and resonance assignments have remained incomplete. Here we present the near-complete (>97 %) backbone assignments of a C145A variant of Mpro (Mpro C145A ) both with, and without, the N-terminal auto-cleavage substrate sequence, in its native homodimeric state. We also present SILLY (Selective Inversion of thioL and Ligand for NOESY), a simple yet effective pseudo-3D NMR experiment that utilizes NOEs to identify interactions between Cys-thiol or aliphatic protons, and their spatially proximate backbone amides in a perdeuterated protein background. High protection against hydrogen exchange is observed for 10 of the 11 thiol groups in Mpro C145A , even those that are partially accessible to solvent. A combination of SILLY methods and high-resolution triple-resonance NMR experiments reveals site-specific interactions between Mpro , its substrate peptides, and other ligands, which present opportunities for competitive binding studies in future drug design efforts.

Real-time Exchange of the Lipid-bound Intermediate and Post-fusion States of the HIV-1 gp41 Ectodomain.

The envelope glycoprotein gp41 of the HIV-1 virus mediates its entry into the host cell. During this process, gp41 undergoes large conformational changes and the energy released in the remodeling events is utilized to overcome the barrier associated with fusing the viral and host membranes. Although the structural intermediates of this fusion process are attractive targets for drug development, no detailed high-resolution structural information or quantitative thermodynamic characterization are available. By measuring the dynamic equilibrium between the lipid-bound intermediate and the post-fusion six-helical bundle (6HB) states of the gp41 ectodomain in the presence of bilayer membrane mimetics, we derived both the reaction kinetics and energies associated with these two states by solution NMR spectroscopy. At equilibrium, an exchange time constant of about 12 seconds at 38 °C is observed, and the post-fusion conformation is energetically more stable than the lipid-bound state by 3.4 kcal mol-1. The temperature dependence of the kinetics indicates that the folding occurs through a high-energy transition state which may resemble a 5HB structure. The energetics and kinetics of gp41 folding in the context of membrane bilayers provide a molecular basis for an improved understanding of viral membrane fusion.

Simultaneous Quantification of H2O2 and Organic Hydroperoxides by 1H NMR Spectroscopy.

Due to similar reactivity of organic hydroperoxides (OHPs), an HPLC separation step is typically required for their indirect (chemical) quantification in mixtures. The high sensitivity of chemical shifts to chemical structure makes NMR an ideal tool for the simultaneous quantification of OHPs in mixtures, but the concentration of these analytes in the samples of interest is usually well below the sensitivity of standard NMR experiments. This sensitivity problem can be mitigated by taking advantage of the fact that the z magnetization of the H2O2 resonance recovers at the rate of hydrogen exchange with water, which is significantly faster than longitudinal relaxation, thus enabling very fast scanning for signal-to-noise enhancement. An adaptation of the E-BURP2 pulse is described that suppresses the water signal by more than 4 orders of magnitude, yielding uniform excitation of peroxide signals without interference of the ca. 108-fold stronger H2O resonance. We demonstrate the method for a mixture of OHPs and report the chemical shifts for multiple OHPs that are of interest in atmospheric chemistry. As shown for hydroxymethyl hydroperoxide, the chemical decay of OHPs can be tracked directly by NMR spectroscopy.

Advances in NMR Spectroscopy of Weakly Aligned Biomolecular Systems.

The measurement and application of residual dipolar couplings (RDCs) in solution NMR studies of biological macromolecules has become well established over the past quarter of a century. Numerous methods for generating the requisite anisotropic orientational molecular distribution have been demonstrated, each with its specific strengths and weaknesses. In parallel, an enormous number of pulse schemes have been introduced to measure the many different types of RDCs, ranging from the most widely measured backbone amide 15N-1H RDCs, to 1H-1H RDCs and couplings between low-γ nuclei. Applications of RDCs range from structure validation and refinement to the determination of relative domain orientations, the measurement of backbone and domain motions, and de novo structure determination. Nevertheless, it appears that the power of the RDC methodology remains underutilized. This review aims to highlight the practical aspects of sample preparation and RDC measurement while describing some of the most straightforward applications that take advantage of the exceptionally precise information contained in such data. Some emphasis will be placed on more recent developments that enable the accurate measurement of RDCs in larger systems, which is key to the ongoing shift in focus of biological NMR spectroscopy from structure determination toward gaining improved understanding of how molecular flexibility drives protein function.

Concordance of X-ray and AlphaFold2 Models of SARS-CoV-2 Main Protease with Residual Dipolar Couplings Measured in Solution.

The 68-kDa homodimeric 3C-like protease of SARS-CoV-2, Mpro (3CLpro/Nsp5), is a promising antiviral drug target. We evaluate the concordance of models generated by the newly introduced AlphaFold2 structure prediction program with residual dipolar couplings (RDCs) measured in solution for 15N-1HN and 13C'-1HN atom pairs. The latter were measured using a new, highly precise TROSY-AntiTROSY Encoded RDC (TATER) experiment. Three sets of AlphaFold2 models were evaluated: (1) MproAF, generated using the standard AlphaFold2 input structural database; (2) MproAFD, where the AlphaFold2 implementation was modified to exclude all candidate template X-ray structures deposited after Jan 1, 2020; and (3) MproAFS, which excluded all structures homologous to coronaviral Mpro. Close agreement between all three sets of AlphaFold models and experimental RDC data is found for most of the protein. For residues in well-defined secondary structure, the agreement decreases somewhat upon Amber relaxation. For these regions, MproAF agreement exceeds that of most high-resolution X-ray structures. Residues from domain 2 that comprise elements of both the active site and the homo-dimerization interface fit less well across all structures. These results indicate novel opportunities for combining experimentation with molecular dynamics simulations, where solution RDCs provide highly precise input for QM/MM simulations of substrate binding/reaction trajectories.

Transient lipid-bound states of spike protein heptad repeats provide insights into SARS-CoV-2 membrane fusion.

Entry of SARS-CoV-2 into a host cell is mediated by spike, a class I viral fusion protein responsible for merging the viral and host cell membranes. Recent studies have revealed atomic-resolution models for both the postfusion 6-helix bundle (6HB) and the prefusion state of spike. However, a mechanistic understanding of the molecular basis for the intervening structural transition, important for the design of fusion inhibitors, has remained elusive. Using nuclear magnetic resonance spectroscopy and other biophysical methods, we demonstrate the presence of α-helical, membrane-bound, intermediate states of spike's heptad repeat (HR1 and HR2) domains that are embedded at the lipid-water interface while in a slow dynamic equilibrium with the postfusion 6HB state. These results support a model where the HR domains lower the large energy barrier associated with membrane fusion by destabilizing the host and viral membranes, while 6HB formation actively drives their fusion by forcing physical proximity.

NMR characterization of H2O2 hydrogen exchange.

Quantification of H2O2 concentration in aqueous solutions is of interest in many fields. It usually is based on indirect methods that rely on oxidation reactions that turn on/off fluorescent probes. Such methods can suffer from reaction incompleteness and interfering chemical species. We describe optimization of NMR detection that enables direct quantification of H2O2 down to the nanomolar range. Taking advantage of fast hydrogen exchange (HX) between H2O2 and water permits the use of very short interscan delays, greatly increasing sensitivity. The specific acid-, base- and water-catalyzed HX rates at 2 °C were measured to be 2.1 × 107, 6.1 × 109, and 1.4 × 10-1 M-1s-1, respectively, which result in a minimum HX rate at pH 6.2. Furthermore, the exchange is accelerated by general acid/base catalysis. MES and phosphate buffers catalyze HX strongest in their unprotonated forms. For imidazole, only the unprotonated form catalyzes HX, which contrasts with acetic acid where only the protonated state catalyzes exchange. Inorganic salts such as sodium chloride and azide have negligible effect on HX. We present optimal conditions for accurate measurement of H2O2 concentrations as low as 40 nM in aqueous samples in a few hours.

Hydrating the respiratory tract: An alternative explanation why masks lower severity of COVID-19.

The seasonality of respiratory diseases has been linked, among other factors, to low outdoor absolute humidity and low indoor relative humidity, which increase evaporation of water in the mucosal lining of the respiratory tract. We demonstrate that normal breathing results in an absorption-desorption cycle inside facemasks, in which supersaturated air is absorbed by the mask fibers during expiration, followed by evaporation during inspiration of dry environmental air. For double-layered cotton masks, which have considerable heat capacity, the temperature of inspired air rises above room temperature, and the effective increase in relative humidity can exceed 100%. We propose that the recently reported, disease-attenuating effect of generic facemasks is dominated by the strong humidity increase of inspired air. This elevated humidity promotes mucociliary clearance of pathogens from the lungs, both before and after an infection of the upper respiratory tract has occurred. Effective mucociliary clearance can delay and reduce infection of the lower respiratory tract, thus mitigating disease severity. This mode of action suggests that masks can benefit the wearer even after an infection in the upper respiratory tract has occurred, complementing the traditional function of masks to limit person-to-person disease transmission. This potential therapeutical use should be studied further.

A lowly populated, transient ß-sheet structure in monomeric Aß1-42 identified by multinuclear NMR of chemical denaturation.

Chemical denaturation is a well-established approach for probing the equilibrium between folded and unfolded states of proteins. We demonstrate applicability of this method to the detection of a small population of a transiently folded structural element in a system that is often considered to be intrinsically fully disordered. The 1HN, 15N, 13Cα, and 13C' chemical shifts of Aß1-40 and Aß1-42 peptides and their M35-oxidized variants were monitored as a function of urea concentration and compared to analogous urea titrations of synthetic pentapeptides of homologous sequence. Fitting of the chemical shift titrations yields a 10 ± 1% population for a structured element at the C-terminus of Aß1-42 that folds with a cooperativity of m = 0.06 kcal/mol·M. The fit also yields the chemical shifts of the folded state and, using a database search, for Aß1-42 these shifts identified an antiparallel intramolecular ß-sheet for residues I32-A42, linked by a type I' ß-turn at G37 and G38. The structure is destabilized by oxidation of M35. Paramagnetic relaxation rates and two previously reported weak, medium-range NOE interactions are consistent with this transient ß-sheet. Introduction of the requisite A42C mutation and tagging with MTSL resulted in a small stabilization of this ß-sheet. Chemical shift analysis suggests a C-terminal ß-sheet may be present in Aß1-40 too, but the turn type at G37 is not type I'. The approach to derive Transient Structure from chemical Denaturation by NMR (TSD-NMR), demonstrated here for Aß peptides, provides a sensitive tool for identifying the presence of lowly populated, transiently ordered elements in proteins that are considered to be intrinsically disordered, and permits extraction of structural data for such elements.

Four-dimensional NOE-NOE spectroscopy of SARS-CoV-2 Main Protease to facilitate resonance assignment and structural analysis.

Resonance assignment and structural studies of larger proteins by nuclear magnetic resonance (NMR) can be challenging when exchange broadening, multiple stable conformations, and 1H back-exchange of the fully deuterated chain pose problems. These difficulties arise for the SARS-CoV-2 Main Protease, a homodimer of 2 × 306 residues. We demonstrate that the combination of four-dimensional (4D) TROSY-NOESY-TROSY spectroscopy and 4D NOESY-NOESY-TROSY spectroscopy provides an effective tool for delineating the 1H-1H dipolar relaxation network. In combination with detailed structural information obtained from prior X-ray crystallography work, such data are particularly useful for extending and validating resonance assignments as well as for probing structural features.

NUScon: a community-driven platform for quantitative evaluation of nonuniform sampling in NMR.

Although the concepts of nonuniform sampling (NUS​​​​​​​) and non-Fourier spectral reconstruction in multidimensional NMR began to emerge 4 decades ago , it is only relatively recently that NUS has become more commonplace. Advantages of NUS include the ability to tailor experiments to reduce data collection time and to improve spectral quality, whether through detection of closely spaced peaks (i.e., "resolution") or peaks of weak intensity (i.e., "sensitivity"). Wider adoption of these methods is the result of improvements in computational performance, a growing abundance and flexibility of software, support from NMR spectrometer vendors, and the increased data sampling demands imposed by higher magnetic fields. However, the identification of best practices still remains a significant and unmet challenge. Unlike the discrete Fourier transform, non-Fourier methods used to reconstruct spectra from NUS data are nonlinear, depend on the complexity and nature of the signals, and lack quantitative or formal theory describing their performance. Seemingly subtle algorithmic differences may lead to significant variabilities in spectral qualities and artifacts. A community-based critical assessment of NUS challenge problems has been initiated, called the "Nonuniform Sampling Contest" (NUScon), with the objective of determining best practices for processing and analyzing NUS experiments. We address this objective by constructing challenges from NMR experiments that we inject with synthetic signals, and we process these challenges using workflows submitted by the community. In the initial rounds of NUScon our aim is to establish objective criteria for evaluating the quality of spectral reconstructions. We present here a software package for performing the quantitative analyses, and we present the results from the first two rounds of NUScon. We discuss the challenges that remain and present a roadmap for continued community-driven development with the ultimate aim of providing best practices in this rapidly evolving field. The NUScon software package and all data from evaluating the challenge problems are hosted on the NMRbox platform.

Concentration-Dependent Structural Transition of the HIV-1 gp41 MPER Peptide into α-Helical Trimers.

The membrane proximal external region (MPER) of HIV-1 gp41 contains epitopes for at least four broadly neutralizing antibodies. Depending on solution conditions and construct design, different structures have been reported for this segment. We show that in aqueous solution the MPER fragment (gp160660-674 ) exists in a monomer-trimer equilibrium with an association constant in the micromolar range. Thermodynamic analysis reveals that the association is exothermic, more favorable in D2 O than H2 O, and increases with ionic strength, indicating hydrophobically driven intermolecular interactions. Circular dichroism, 13 Cα chemical shifts, NOE, and hydrogen exchange rates reveal that MPER undergoes a structural transition from predominately unfolded monomer at low concentrations to an α-helical trimer at high concentrations. This result has implications for antibody recognition of MPER prior to and during the process where gp41 switches from a pre-hairpin intermediate to its post-fusion 6-helical bundle state.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy.

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (Rh) of biomolecules in solution. For unfolded proteins, the Rh provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the Rh of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by 13C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that Rh values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the Rh of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Modulating the Stiffness of the Myosin VI Single α-Helical Domain.

Highly charged, single α-helical (SAH) domains contain a high percentage of Arg, Lys, and Glu residues. Their dynamic salt bridge pairing creates the exceptional stiffness of these helical rods, with a persistence length of more than 200 Å for the myosin VI SAH domain. With the aim of modulating the stiffness of the helical structure, we investigated the effect, using NMR spectroscopy, of substituting key charged Arg, Lys, Glu, and Asp residues by Gly or His. Results indicate that such mutations result in the transient breaking of the helix at the site of mutation but with noticeable impact on amide hydrogen exchange rates extending as far as ±2 helical turns, pointing to a substantial degree of cooperativity in SAH stability. Whereas a single Gly substitution caused transient breaks ∼20% of the time, two consecutive Gly substitutions break the helix ∼65% of the time. NMR relaxation measurements indicate that the exchange rate between an intact and a broken helix is fast (>300,000 s-1) and that for the wild-type sequence, the finite persistence length is dominated by thermal fluctuations of backbone torsion angles and H-bond lengths, not by transient helix breaking. The double mutation D27H/E28H causes a pH-dependent fraction of helix disruption, in which the helix breakage increases from 26% at pH 7.5 to 53% at pH 5.5. The ability to modulate helical integrity by pH may enable incorporation of externally tunable dynamic components in the design of molecular machines.